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Abstract 

Atom-metal  surface  overlap  integrals  are  of  utmost  importance  in  surface 
energy  calculations.  Direct  numerical  evaluation  of  these  triple  integrals  can 
be  very  time  consuming.  However,  we  have  developed  an  exact  algebraic 
expression,  where  formulas  for  the  coefficients  are  given  for  both  the  general 
case  and  the  special  case  where  the  parallel  wavevector  is  zero.  Some 


numerical  examples  of  the  overlap  for  H  on  A1  are  given. 
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1.  Introduction 

The  interaction  of  an  adatom  with  a  metal  surface  is  of  utmost  importance  in 
surface  physics  (and  surface  science  in  general).  Several  researchers1  ^  have 
expended  a  great  deal  of  effort  in  its  determination.  To  evaluate  this 
interaction  potential,  one  common  method  is  to  first  expand  the  total 
wavefunction  in  a  mixed  basis  set, 

♦  ,(r)  -  l  c*  $a(r)  +  I  cj<$k^  *  O) 

a  k 

where  4>a(r)  are  wavefunctions  that  are  localized  on  the  adatom  with  quantum 

number  a,  and  4>k(r)  are  the  delocalized  wavefunctions  of  the  metal  with  quantum 

number  k.  One  of  the  difficulties  involved  with  such  an  expansion  is  the  need 

14 

to  remove  the  overcompleteness.  Lundqvist  has  suggested  that  this  can  best  be 
done  by  the  requirement 

I  \  Sa.k  '  0  •  (2> 

where  the  overlap  integral  is  given  by 

sa,k  "  <alk>  *  (3 ) 

Consequently,  a  knowledge  of  this  overlap  is  needed  to  remove  the 
overcompleteness.  Furthermore,  to  solve  the  secular  determinant  for  Eq.  (1), 
one  needs  to  know  the  values  of  the  overlap  integral  and  two  interaction 
Integrals.  The  interaction  integrals  in  turn  can  be  related  to  each  other  and 
the  overlap  via  the  Hermitian  property 

Ea  \,a  *  ‘  Ek  Sa.k  *  <alValk>  •  <*> 


where  E  and  V  are  the  eigenvalue  and  potential  of  the  isolated  atom,  and  E 
and  are  the  eigenvalue  and  potential  of  the  isolated  metal.  Consequently,  an 
evaluation  of  the  overlap  is  vital  for  solving  the  secular  deteminant. 

In  this  paper,  we  show  that  a  closed  form  expression  for  the  overlap  can  be 
obtained.  We  define  our  system  in  the  next  section,  and  following  this  we  give 
the  details  of  evaluating  S  .  Finally,  v  „>  discuss  our  results  and  suggest 

3  i  K 

further  uses  of  this  calculation. 

2.  The  System 

We  shall  examine  the  case  of  a  single  atom  impinging  on  a  simple  metal 
surface.  The  electrons  associated  with  the  atom  will  be  approximated  by  Slater 

7 

orbitals, 


$  ( r)  « 
a 


[Uilf— ] 

1  (2n) !  J 


1/2 


n-1  -£r  vm, .  . 

r  e  Ya(e,<}>) 


where  Y^(e,<t)  are  the  spherical  harmonics,  a*(n,£,m)  are 
numbers,  and  E,  is  the  orbital  factor.  The  origin  of  our 
on  the  surface;  therefore. 


(5) 

the  atomic  quantum 
system  is  chosen  to  be 


r  ■=  v  r^  +  (z-z  )2  ,  (6) 

where  z  and  r((  are  the  position  of  the  electron  perpendicular  and  parallel  to 
the  surface,  and  z  is  the  distance  of  the  atom  from  the  surface.  The  Slater 

d 

orbital  choice  was  made  for  computational  efficiency.  Furthermore,  if  V  is 

a 

Coulombic,  the  interaction  integral  for  the  atomic  potential  will  be 


.  2  a _ s  . 

v/(2n)  (2n-1 )  a',k 


(7) 


where  Z  is  the  nuclear  charge  and  a'  -  (n-1,£,m).  Consequently,  all  integrals 


of  interest  can  be  written  in  terms  of  the  overlap. 


5?.  V. 


Previously,  the  metal  has  been  modeled  within  the  truncated  jellium 
1-5 

approximation.  Instead  of  using  the  numerical  wavefunctions  from  this 

model,  we  shall  consider  the  metal  electrons  as  particle-in-a-box.  Such 

wavefunctions  are  good  approximations  to  the  jellium  model  and  provide  the 

basis  set  necessary  for  a  more  exact  approach  within  the  nearly-free-electron 
8 

approximation.  These  particles-in-a-box  wavefunctions  are 

■+  ■* 

■*  -1/2  *r 

♦,  (r)  *  A  1  e  f ( z)  ,  (8) 

k  s 

where  A  is  the  surface  normalization  area,  and  k  is  the  component  of  the 

S  11 

electronic  wavevector  parallel  to  the  surface.  f(z)  is  the  one-dimensional 
particle-in-a-box  wavefunction: 


f(z)  -  (jO  sin(kzz  +  6k) 


for  -L  <  z  <  0, 


f(z)  -  (£)  sinek  e_qz 


for  z  >  0,  and 


f(z)  -  (£)  sin(-kzL  ♦  ek)  eq(z+L) 


for  z  <  -L,  with 


2W  -  k 


where  L  is  the  thickness  of  the  metal,  kz  is  the  component  of  the  electronic  wave- 
vector  perpendicular  to  the  surface,  and  W  is  the  sum  of  the  work  function  and 
the  Fermi  energy. 

To  calculate  the  overlap  of  the  wavefunctions  given  in  Eqs.  (5)  and  (8),  we 
must  solve  the  integral 


1/2 


5 


<,  .-1/2  r(2U2nM 

‘"a.k  s  1  <2n) ! 


n-1  -'r  vm*„  . 

d  r  r  e  o ,  <? )  e 


ik  *r 


f(z+z  )  ,  (12) 

o 


where  the  zg  is  needed  in  the  metal  wavefunction  since  our  origin  is  centered  on 
the  adatom.  At  first  this  integral  appears  quite  complicated,  due  to  the  broken 
symmetry  caused  by  the  surface.  However,  we  can  carry  out  an  evaluation  by 
means  of  a  Fourier  expansion. 

3.  Evaluation 

The  metal  wavefunction,  Eq.  (8),  can  be  represented  by  a  Fourier  transform, 


*k(r)  «=  ds  ik(s)  eis’r 


where  s  is  the  Fourier  coordinate.  The  transform  of  the  wavefunction  is  given 


m 


f(s  ) 

Va)  -  -7 h 

a  #> 


.  IS  z 

f(sz)  =  —  dz  e  z  f(z+za)  , 


where  6(k(i-s|()  is  the  Dirac  delta  function.  Using  Eq.  (9)  in  Eq.  (15),  we  obtain 

,  _  i /o  *  is  (z-z  ) 

f(s  )  =  — — (-)  { sin6  j  dz  e  Z  a  e  qZ 

z  2r  L  k  > 


-L  is  (z-z  ) 


+  sin  (S.  -k  L)  \  dz  e 
k  z  / 


a  eq(Z+L)  +  J  dz  e 


0  is  (z-z  ) 


sin(k  z+6  ) :  (16, 
z  k 


Evaluating  these  integrals  is  straightforward  and  gives 

-is  z  -is  (z  +L) 

.  -1/2  i  sine,  e  i  sin( 6. -k  L)  e  z  a 

<?>  i— rVi^ - - 

z  M  z  w 


19k  -lek  l[(V.t)L.9k] 

[  £ —  ♦  £ - - - £ - —  1 1 

1  k  *s  k  -s  k  +s  k  -s  J 1 

2  Z  Z  Z 


k  +s 
z  z 


k  -s 
z  z 


V  •  ■  V  V  v  h *■  .  •  .  *  .  *  *  - . 
1  'mV*/'  V 


Combining  Eqs.  (12)  and  (13),  we  can  express  the  overlap  as 


o 


a,k 


n-1  -£r  , .  1st 

dr  r  e  (0.4»)  e 


(18  ) 


The  exponential  can  be  expanded  in  terms  of  the  spherical  harmonics, 
•  V 


■¥  + 
1st 


*•»  I  I  i  C  (es.^)Y™,(e,0)j,f(sr)  , 

i'-O  m*—  V  S 


(19) 


where  j£,(sr)  is  the  spherical  Bessel  function.  Using  this  in  Eq.  (IE)  and  the 
orthonormality  of  the  spherical  harmonics,  the  overlap  reduces  to 


a,k 


(2E)20*1  1 

.  (2n) !  J 


1  /2 


•4tt  i 


ds 


4>(  s) 


<vv 


dr  r 


n+1 


j£(sr) 


(20 


By  Eq.  CU),  this  can  be  further  reduced  to 

1/2 

2n*1 


a,k  *■  4ir  ( X.+ 1  in  j  ) !  (2n) !  ■*  s 


-1/2 


'irad1,. 


*  e 


asz  f<sz>  PIn,|<SS) 


dr  rn+1  j£(sr)  e  , 


(21) 


where  <tk  is  the  azimuthal  angle  in  k-space, 


is  the  associated  Legendre 


function,  and  s  is  now  given  by 


,/~2  2 

V  k  +  s 


(22) 


The  spherical  Bessel  function  can  readily  be  expanded  as 


10 


p  t — 1 

1  r  1  ((H)!  r  isr 

t(Sr)  '  2s?  -  t  ^ 

t-0  t!(£-t)!(2sr)C 


+  e-isr(_1)t-£-1 


(23) 


Using  this  expression,  we  can  evaluate  the  inner  integral  in  Eq.  (21): 


c  „-1/2 

k  n  T  1 

( i ♦ t ) ! ( n- 1 ) ! 

s  *  A  e 

a,k  s 

D  l 

t-0 

t!  ( ft-t ) ! 

«  dsz  f(sz)  p[ 

. 

1  S 

1 

8  (2s)t+1 

If  ’ 

(5-is)n-t4 

where  for  convenience, 

we  define 

1/2 

„  .  ■..ft?**')  (HH)-'  UP  1 

j-0  Cn-t+i-J>!  (C2^t^)n-t+1 


Eq.  (2fc)  can  be  solved  by  considering  the  integral  as  extended  over  a  semicircle 

in  the  negative  complex  plane.  The  total  integral  over  this  enclosed  path  gives 
+«  . 

(  ds  ...  +  ids  ...  =  -2iri  7  R.  ,  (27) 

_oc  c  g  J 

where  the  integrand  is  the  same  as  in  Eq.  (26).  The  C  in  the  second  integral 

implies  integration  over  the  semicircle  part  of  the  curve,  and  R_.  are  the  residues 
of  the  poles  contained  within  the  total  curve.  Since  the  integral  is  in  the 
negative  complex  plane,  the  exponentials  in  f  (s^)  goes  to  zero  at  infinity; 
the  second  integral  vanishes.  Consequently,  the  integral  in  Eq.  (26)  only 
depends  on  the  poles  in  the  negative  complex  plane.  Inspection  of 
the  integrand  reveals  three  such  poles:  (1)  a  pole  of  order  1  at  s  -  -iq 

/ ~2  5 

contained  in  f(s  );  (2)  a  pole  of  order  n+1  at  s  -  -We,  +k..  from  the  fraction 

Z  Z  II 

in  the  sum;  and  (3)  a  pole  of  order  t  at  s  ■  -ik(|  due  partially  to  the  Legendre 
function  and  partically  to  the  fraction  it  multiplies. 

Using  Eq.  (17),  we  can  easily  determine  the  residue  at  the  first  pole: 


where  we  have  ignored  terms  that  vanish  as  the  metal  thickness,  L,  becomes 


large.  In  order  to  evaluate  the  other  two  residues,  we  first  define 


i  sine  *9k  i6k 

F!ls,>  ■  [rrrr  -  ;  <rnr  4  rrr>l 


s  +  ia.  2  k  +s  k  -s 

Z  K  Z  Z  Z  Z 


£  ,t-1 


I  1  (£>t)!  (r.-t)!  |m|  fz.  r 

t-0  t!  ( £-t) !  £  vs  ‘  L 


(2s) 


t+1 


]  [ 


— i - .  ( 

(C-is)^1  (CMs)n~t+1 


n  2 

Consequently,  we  can  now  write  the  residue  at  5  «  -i/£  +  k 


as 


r  .  —  t -2_) 
2  2it  'LA  ’ 
3 


1/2  ~i  mt 


K  »_  2l|[s  „^Jin4'  FUs  ,  ,Val 


t2.k2-  < 

t  K(| 


If  we  let  x  -  is  ,  then 
z 


n+1 


-  k  <£>  k  {l-'Rl  F<*>  H  „ 

s  dx  t  )  ,/.2^,2 

J  +k„ 


x-vC 


II 
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Using  Leibnitz'  formula  for  the  differentiation  of  the  product,  we  obtain 


2  2ir  vLA 

s 


n!  h»0  (n-h>!  h!  a 


n+1 


F(x) 

dx  1 


v42+k2 


i 


tf  we  now  define 

n-h+1 „  .h 
d 

dx 


B,  = 


(-1) 


h  (n-h) 


_ ’  si  \  U, 

'■  •>!  „vh  ( 


n+1 


2+k»  F(x)l  m 

11  3  x-vr -k; 


we  can  then  write  the  second  residue  as 


1,2,  ‘lm*k  .  5  -  n-h,  za 

r,  (CT>  e  <J  „Va  >  e 

s  h=0 


Similarly,  we  can  obtain  the  residue  at  s  -  -ik  as 

z  II 


1-1 


i  ,  2  , 
r:  (,-T")  e 


.  ,  ,  p  „  1-1  ~hv  kllza 
imc^  (  2  Chzg  )  e 


3  2  it  LA 


h=0 


(32 


(3£: 


(34: 


(35 


for  1  >  0.  If  1  -  0,  there  is  no  residue  at  s  -  -ik  .  The  coefficients  are 

z  II 


given  by 


-  .  <-ui_S._  si  irx.k  ji  F(s)| 

'h  (1-l-h)!  h!  dxh  }lx  llJ  u;j 


(36! 


x«k 


II 


In  the  spirit  of  the  second  and  third  residues,  we  can  rewrite  the  first,  Eq. 
(26),  as 


p  'Im*k  .  /QZa 

R1  “  2tt  (LA  )  e  A 


(37 


where 


A  ■  D  sin6  V  ?\m\(  ~q— \ - J - 


Izll 


t-JL-1 


g  *  g  - 


n-t+1 


(36) 


Since  we  now  have  the  value  of  the  residue  at  each  pole,  we  are  able  to  obtain 
the  overlap 


3 

S.  k  *  -2*i  I  Rt  . 

a,k  .=1  i 


(39 ) 


where  the  minus  sign  is  due  to  our  integrating  in  the  clockwise  direction. 
Using  Eqs.  (3H),  (35)  and  (37).  we  obtain 


a.k 


’'2  -InO  -qz  r, 

’lF1  e  l#e  4  lJnBh  za  1 

s  h=0 


£  4kll  za 


rV  l-l -hi  kllzai 
l  ' ln  Ch  za  1  e  1 
n*0 


(^0) 


Thus,  we  have  reduced  the  overlap  integral  to  a  simple  algebraic  expression  in 
terms  of  the  distance  of  the  adatom  from  the  surface.  Furthermore,  as  we  shall 
see  below,  A,  and  are  real;  consequently,  all  imaginary  contributions  to 
the  overlap  are  determined  by  the  simple  phase  factor  in  Eq.  (HO). 

<< .  Numeric  •'edures 

In  order  to  c  oute  values  for  the  overlap,  Eq.  (MO),  one  must  solve  for  the 
coefficients.  Th  procedure  for  A  is  straightforward,  and  its  value  is  given  by 
Eq.  (3£i).  ana  C^,  however,  dependent  on  several  differentiations  of  the 


function  F(is  ).  One  simplification  is  obtained  by  observing  that  F(is  )  is 
z  z 

completely  real.  This  can  be  shown  by  making  the  substitution  x  *  is  in  Eq. 

z> 

(29): 


F(x)  -  [■ 


k  cose,  -  x  sine, 

z _ k _ k_ 

2^,  2 
x  +k 

z 


1  y  U+t)!  (n-t)|  rl m| /  -_x_ 

J0  t,  (i-t),  I 


.  t-£-l 


[$>/?- k^)  (i +  j x2_k(i) 


n-t+l 


fc  -  /x2-kjj. 


n-t+l 


To  proceed  with  the  evaluation  of  B.  and  C.  ,  one  could  obtain  a  closed-form 

h  h 

expression  by  repeated  use  of  Leibnitz’  formula  or  attempt  to  use  a  numerical 
approach  based  on  finite  differences.^  However,  a  simplified  form  of  the  overlap 
can  be  obtained  for  the  special  case  of  k)(  =  0.  Since  the  overlap  would  be 
expected  to  be  greatest  for  large  (k(1  small),  one  could  use  this  special  case 
as  a  basis  for  obtaining  a  general  solution. 

If  we  assume  that  k()  =  0,  the  associated  Legendre  function  in  Eqs.  (3fj)  and 
(HI)  will  vanish  unless  m  =  0.  Consequently,  S  -  0  for  k  -  0  and  m  »  0. 

3  9  K  II 

For  the  case  where  m  =  0,  to  evaluate  one  must  include  in  Eq.  (29)  the  terms 

that  depend  on  L  from  Eq.  (17)  since  these  terms  do  not  vanish  at  s  =0. 

z 

However,  the  value  of  these  terms  assures  that  F(is  )  and,  consequently,  C  goes 

z  n 

to  zero.  Therefore,  for  k.  =  0  and  m  -  0,  the  overlap  becomes 


1/2  ~qz  n 

Sa.k  -‘ET’  I*  e  *  I  £ 
s  h«0 


f  'i  \  Ch]  • 


where  A  can  now  be  simplified  to  the  form 


D  sine  \  ami  Suzlll  JzuL  [ _ 1 _ 

k  ^  t.  (t-t)l  (2q)t+1  lu  +  q)n-t+1 


IzllUli 


U~q) 


n-t+l  J 


Furthermore,  using  Eqs.  (3S)  and  (HI),  reduces  to 


K 


ir 


.  £-h 


B  .  ±1111.  D  ?  iimiinzll 

h  (n-h)!  h!  . ^-t!  (l-t)! 


.h  r  sine. 

4-  |( - iL  + 

.  h  l  x-q 
dx  ^ 


t-0 


k2cos?k  -  x  sinek  (£.x)t  i 

x2+k2  (2x)t+lJx“^  ‘ 


(44) 


Another  application  of  Leibnitz'  formula  leads  to 


Bh-  LZ 


,l-h  „  U.h) 


1)-  _  D  y+tj!  (n-t)!  (-1) 

(n-h)!  tfQ  t!  ( t-t) !  (h-t)! 


.h-t  r  sin0,  k  cose,  -  x  sine, 
d _ [  ( _ k  +  z  k _ \  _ 

„  h-tLv  x-q  2  2  ‘  ._ 

dx  M  x  +k  (2x 


hr]  • 

)  x-C 


(4  5) 


where  (!,h)  implies  using  either  £,  or  h,  whichever  is  smaller. 

We  have  evaluated  the  overlap  integral,  Eq.  (42),  for  a  variety  of 
interacting  orbitals  for  hydrogen  atom  on  the  surface  of  aluminum,  and  the 
results  are  depicted  in  Fig.  1.  As  can  be  seen  from  this  figure,  the  overlap 
becomes  large  as  n  becomes  large,  which  is  due  to  the  fact  that  the  overlap  does 
not  become  a  maximum  until  the  volume  occupied  by  the  electronic  wavefunction  of 
the  atom  approaches  the  size  of  the  metal.  This  happens  at  the  point  of 
dissociation,  n  ■*  <*>. 

The  peculiar  structure  of  the  2s  overlap  is  caused  because  q  *  £,  i.e.,  at 
the  exact  point  of  equality  A  cancels  bQ.  Consequently,  if  the  atom  is  on  the 
surface  (z  «  0),  the  overlap  is  zero.  For  q  <  5,  the  metal  surface 
damping  dominates  the  overlap  at  small  z;  for  q  >  the  atomic  damping 


dominates. 
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5.  Discussion 


We  have  demonstrated  that  the  overlap  Integrals  for  the  wavef unctions  of  an 
atom  with  a  metal  surface  can  be  written  as  algebraic  expressions  In  term  of  the 
distance  from  the  surface.  Furthermore,  we  have  shown  that  these  expressions 
are  real  except  for  a  simple  multiplicative  phase  factor.  Analytical 
expressions  for  the  coefficients  of  the  overlap  formulas  have  also  been  pre¬ 
sented.  Solving  these  closed-form  expressions  for  the  overlaps  will  be  much 
more  computationally  efficient  than  direct  numerical  solution  of  the  overlap 
integrals.  Since  other  integrals  of  interest  in  surface  calculations  can  be 
written  in  terms  of  the  overlaps,  this  efficiency  will  be  transferred  to  atom- 
metal  surface  potential  energy  calculations.  Such  integrals  are  also  extremely 
useful  in  studying  such  problems  as  atom-surface  charge  transfer,  scattering  and 
energy  coupling.  These  problems  are  of  ongoing  interest  in  our  research. 
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1.  The  overlap  Sa>2(z  )  in  units  of  the  normalization  factor,  (2/LA  ) 

versus  the  dlstancl  from  the  surface,  z  ,  in  atomic  units.  The  data  are  To 
hydrogen  s  orbitals  overlapping  the  FerSi  wavefunction  of  aluminum.  The 
curves  are  labelled  by  the  principal  quantum  numbers. 
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